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Abstract 

Ab initio studies of structural, elastic and electronic properties of the tetrag- 
onal perovskite-type PbZro.5Tio.5O3 are presented using the pseudo-potential 
plane wave method within the density functional theory in generalized gra- 
dient approximation. The calculated equilibrium lattice parameters remain 
in a good agreement with the available experimental data. The bulk mod- 
ulus obtained from the Birch-Murnaghan equation of state is calculated as 
Bo=189 GPa, and the gap energy Eg=2.1 eV. The largest differences between 
calculated and nominal charges exist for Zr atoms (the biggest value is about 
50% larger than +4). The influence of the strain in the 0.1%-0.3% range was 
studied. 
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1. Introduction 

Nowadays the hetero-modulus ceramic composites with different mechani- 
cal properties are one of the main directions of the development of new tech- 
nical materials. The most common material is the lead zirconate titanate 
PbZri_2^Ti^.03 (PZT) family which is well-known to have high piezoelectric 
strain constant, electric permittivity and electromechanical coupling It 



Preprint submitted to Computational Materials Science 



October 17, 2012 



is too early to generalize conclusions about using PZT in technical applica- 
tions. They are very wide, starting from crystals, across ceramics, fibres, 
ending on nanostructures and nanoparticles. PZT thin fibres have recently 
attracted considerable attention due to their applications in precision Mi- 
croelectromechanical Systems (MEMS) or integrated electro-ceramic MEMS 
platforms They have a large potential in sensoric, actuatoric and ultra- 
sonic transducer applications [3|. Fibres PZT have a potential for utilization 
in high performance hydrophones and ultrasonic transducer applications be- 
cause the fine fibrous geometry offers the possibility of composite fabrication, 
where damping and reinforcement is combined 4]- PZT fibres are used to 
form active elements within a functional devices as the active fibre composite 
(AFC) [51. 

For better use of PZT advantages the further theoretical and experimental 
investigations are highly required. Ab initio calculations offer one of the most 
powerful tools for carrying out theoretical studies of these properties. First 
ab initio calculations for PbZro.5Tio.5O3 were performed by G. Saghi-Szabo et 
al. 6| . Apart from the stable ferroelectric ground state they determined the 
bulk spontaneous polarization, dynamical charges, and piezoelectric stress 
tensor elements. Piezoelectric parameters achieve maximal values near the 
morphotropic boundary. Experimental results show that up to x = 0.5 of 
the Zr content these fibres maximally contain about 90% the volume fraction 
of a tetragonal phase and 10% a rhombohedral phase in a volume content of 
PZT [7]. So for this reason in this paper we have analyzed structural, elastic, 
electronic properties of the tetragonal phase of PbZro.5Tio.5O3 from ab initio 
calculations, near the morphotropic boundary i. e. close to a; = 0.53 
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Figure 1: The model unit cell with the positions of atoms in PbZro.5Tio.5O3 

The paper is organized as follows: the method of calculations is described 
in Section 2. In Section 3 the results and discussions of the structural, elastic 
and electronic properties as the results of our calculations are presented. 
Finally, the conclusions and remarks are given in Section HI 



2. Method of calculations 



Density functional theory (DFT) 8| calculations were performed in the 
generalized-gradient approximation (GGA) [9]. The localized atomic-like or- 
bitals for the basis set expansion were used as implemented in the SIESTA 



software package 
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Afterinitially carried out tests for different possi- 



bilities the plane wave energy cutoff for the basis set at 400 Ry and a 4 x 4 x 4 
Monkhorst-Pack grid for the Brillouin zone integration were chosen jl2|. The 
unit cell containing 40 atoms is shown in Figure [TJ 

In order to achieve better performance, the most non-valence electrons 



Table 1: Values of calculated and experimental lattice parameters of PbZro.5Tio.5O3 
Lattice parameters Present Work Experiment [14] 

a (A) 4.101 4.042 

c (A) 4.091 4.127 

were replaced by effective norm-conserving Troullier-Martins pseudopoten- 



tials in the Perdew-Burke-Ernzerhof (PBE) parametrization [13|]. In con- 
trast, the electrons Ti of 3d^ is^, Zr id^ 5s^, Pb 4/ Sfi^^ Gs^ 6p^ and O 
2s^ 2p'^ were treated as valence electrons. 

The geometric optimization of atomic positions and lattice parameters 
with respect to the ground state energy was performed using a conjugate- 
gradients method. The obtained values of lattice parameters, in the relax- 
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ation procedure, are very close to the experimental ones 
presented in Table [TJ 

3. Results of calculations 

3.1. Elastic properties 

We have calculated elastic properties both the uniform volume compress- 
ibility (the bulk modulus) as well as elastic constants, which describe different 
directional deformations. The bulk modulus was calculated by fitting the de- 
pendence of the total energy of the unit cell volume to a Birch- Murnaghan 
equation of state [l5| (see Fig. [2]). The dependence of the energy on volume 
was obtained by series of calculations, where the unit cell volume was varied 
between -1 % to +1 %. Such relative volume change is equivalent to the 
hydrostatic pressure in the range from -6 GPa to +6 GPa. From Fig. |5]one 
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Figure 2: The dependence of the total energy on the unit cell volume of PbZro.5Tio.5O3 

can see that the minimum of the total energy Eq = - 22741 eV occurs for 
Vo=549 A'^. The calculated bulk modulus amounts to Bq =189 GPa. 

The dependence of an, 022^ ^"33, (744, 0-55, age, o"i2, a23, and 0-13 on the strain 
was calculated in the range 0.03%-0.3% by applying small strains in selected 
directions with values of the corresponding stress tensor components from 
the Kohn-Sham total energy [8]. The stresses an, O22 and a44, ass are equal 
in pairs due to the tetragonal symmetry in PbZr0.sTi0.sO3. Results are pre- 
sented in Figs[3t^a) and[3t^b). 

Values of an and 022 are about 15% higher than a33. It means that in 
the (a, 6) plane the tetragonal PbZr0.sTi0.sO3 is a little harder than in the c 
direction. The stress a\2 and ai3 describe properties of the PbZr0.sTi0.sO3 in 
the (a, 6) and (a, c) plane, respectively. There are large differences between 
them and other elastic constants. The stress a44 corresponds to the shear de- 
formation in the (6, c) plane. All elastic constants for the shear deformation 
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Figure 3: Calculated directional stress-strain dependences of cth, (T22, cas (a) and (b) (744, 
CT55, CT66, (712 , (J23, and (713 for PbZro.5Tio.5O3 

are significantly smaller than the elastic constants for the uniaxial deforma- 
tion. During calculations we have observed almost the same behavior for 
the compression and the elongation in the dependence of these elastic con- 
stants on the applied strain. Similarly to an experiment 16| one can observe 
the nonlinear behaviors for small strains near 0.05% and next near 0.20%. 
The value 0.25% of the strain is critical - about this value real crystals are 
destroyed. This fact is in very good agreement with the experimental one 
[igI . Basing on these results for the strain range of 0.03%-0.3% we have 
calculated the respective elastic constants, the directional Young's modulus 
E and average values of Young's modulus In FigHlthe dependence of 
the average Young's modulus E on the strain for PbZro.5Tio.5O3 is shown. 
Between strain values of 0.07% and 0.15% there exists only very short region, 
where the Hooke law is valid. 



3.2. Electronic structure 

The calculated total density of states for the tetragonal phase PbZro.5Tio.5O3 
is presented in Figure O The value of a band gap amounts to E^ = 2.1 eV. 
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Figure 4: Calculated values for the average Young's modulus as the function of the strain 
in PbZro.5Tio.5O3 

The net atomic charges, collected in Table [21 indicate the substantial ionic- 
ity of the atoms. The net charges were calculated by the integration of the 
charge density over a sphere of the selected radius. They are the real charges 
in the nodes of the lattice. The integration radii were taken as 1.25 a.u., 1.5 
a.u., 1.7 a.u., and 2.5 a.u. for Ti, Zr, Pb and O atoms, respectively. These 
radia are about a half of the bond length between respective atoms. The 
largest differences between the calculated and nominal charges exist for Zr 
atoms, being about 50% larger than +4. The differences between calculated 
and nominal charges existing for all atoms are supposed to be due to a sub- 
stantial hybridization between the O 2p orbital with valence orbitals of other 
ions. The origin of so large charge differences and physical consequences will 
be further studied. 
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Figure 5: Total density of states for the tetragonal phase PbZro.5Tio.5O3 
Table 2: Calculated net ionic charges in PbZro.5Tio.5O3 







Zr^^ 




Q-II 


Radius of integration (a.u.) 


1.7 


1.5 


1.25 


2.5 


Calculated atomic charge 


3.25 


6.17 


5.77 


-1.24 


Ionic radius (pm) 


119^^ 77.5^^ 


72^^ 


67^^^ 60.5-^^ 


140"^^ 



4. Discussion and conclusions 

From the analysis of the stress-strain directional dependences and the 
Young modulus-strain dependency one can see that in PbZro.5Tio.5O3 near 
morphotropic boundary the Hooke law is carried out only on the range of 
the strain 0.07%-0.15%. Values of directional Young's modulus are approx- 
imately constant in this range. Such behaviors seem to be typical for the 
stress-strain behavior of PbZro.5Tio.5O3 fibres as has been written in Refs 
jl^ [igI . [i8 |. Theoretical results are over estimate in comparison with the ex- 
perimental ones - it is supposed to be due to defects occurring even in the best 
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real crystals (calculations have been performed for the stoichiometric struc- 
ture and the temperature at K). The experimental value of the Young's 



modulus along the radius 
in the room temperature 



'or ceramic PZT fibres are about 108 GPa-120 GPa 
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The same effect 
perovskite crystal 



las been observed in ab initio calculations for other 
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2l|. Knowing from theoretical calculations that point 
defects cause decreasing almost by a half of calculated elastic constants of 
the ideal crystal it is necessary to compare the calculated values of these 
parameters for the defected crystals with the experimental values of elastic 
constants for tetragonal PZT measured at low temperatures. 
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